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Abstract 

Network motifs are characteristic patterns which occur in the networks essentially more 
frequently than the other patterns. For five motifs found in S. Itzkovitz, U. Alon, Phys. 
Rev. E, 2005, 71, 026117-1, hierarchical random graph models are proposed in which the 
motifs appear at each hierarchical level. A rigorous construction of such graphs is performed 
and a number of their structural properties are analyzed. This includes degree distribution, 
amenability, clustering, and the small world property. For one of the motifs, annealed phase 
transitions in the Ising model based on the corresponding graph are also studied. 


1 Introduction 

In view of the complexity and unknown organizing principles of large real-world networks, they 
usually are modeled by means of random graphs, the study of which traces back to P. Erdos 
and A. Renyi [1]. Many of such networks contain characteristic patterns recurring essentially 
more frequently than the other ones. These are network motifs ISIEISIIS]. Quite often real 
networks are build up mostly of motifs, which thus can be treated as constructing units for 
their modeling, cf. [5]. In [3], the authors introduced a random graph model based on some 
geometric principles (constraints). Then they compared the appearance of eight elementary 
three- and four-node patterns in their model with the same characteristics of the Erdos-Renyi 
random graph. It turned out that five of these patterns are motifs for their model, but not for 
the Erdos-Renyi random graph, see Fig. [TJ 



Figure 1: Three and four node motifs Mi, M 2 , M^, M 4 , found in [3]. 

One of the ways to get information about infinite graphs is to study the properties of certain 
models of statistical physics defined thereon. The most popular ones are the Ising and Potts 
models, see [B]. On the other hand, graphs are employed to mimic crystal lattices. For graphs 
with special structure, the critical behavior of the Ising model can be described in an explicit 
and rigorous way. This, in particular, holds for the so-called hierarchical lattices introduced in 
mu- Such lattices are constructed in an algorithmic way by means of basic patterns, e.g., by 
a ‘diamond’, see M 3 in Fig. [TJ The relative simplicity of the theory makes hierarchical lattices 
attractive in studying critical point behavior of various types, see quite recent works laiio] and 
the references therein. A mathematical description of the Gibbs states of the Ising model on such 
graphs was done by P.M. Bleher and E. Zalys in mm- M. Hinczewski and A. Nihat Berker 
m studied the critical point properties of the Ising model on the diamond hierarchical lattice 

‘decorated’ by adding random bonds. In the present paper, we follow the way suggested in [13] 

and introduce hierarchical graphs constructed by means of the motifs shown in Fig. (TJ decorated 
by random bonds which somehow repeat the corresponding motif. We analyze a number of 
their characteristics, such as the average degree, the node degree distribution, amenability, the 
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small-world property. We also study a ferromagnetic phase transition in the Ising model dehned 
on the graph based on Mi. A preliminary study of the models introduced here was performed 

in [UllIS]. 

2 The Graphs: Construction and Structural Properties 

2.1 The construction: informal description 

As is typical for hierarchical graphs, e.g., for hierarchical lattices in [Ellis], the construction is 
carried out in an algorithmic way: at A:-th level, /c G N, one produces a subgraph, say A^,, which 
is then used as a construction element for producing A^+i. The procedure is the same at each 
level. The starting element at level 1 is obtained from the corresponding motif. Let us illustrate 
this in the simplest case based on Mi - the triangle. To obtain Ai, we label the nodes of Mi 
by a, 6, and c, as shown in Fig. [2l The graph A 2 is created in two step. First we take three 
graphs of level 1 and label them by A“, A^ and A^. Thereafter, the triangles are being glued up 
according to the following rule: for i,j G {a, 6, c}, i / j, node i of triangle A{ is glued up with 
node j of triangle A^ The nodes i od triangle A^ remain untouched. These are the external 
nodes of A 2 . The remaining nodes are called internal. The bonds of the initial triangles A^, 
i G {a,b,c} turn into the bonds of A 2 . We call them basic bonds; they are depicted as solid 
lines. At the second step, we add bonds connecting the external nodes in the same way as it is 
in the motif Mi. Such bonds are depicted as dotted lines and called decorations. As a result, 
we obtain the graph A 2 , which has nine basic bonds and three decorations, three external and 
three internal nodes. 


c 



Figure 2: Construction of the graph A 2 based on Mi 


To obtain A^, k = 3,4,..., we repeat the same procedure - take three copies of Ak-i and 
label them by A^_j^, A^_j^, and A|_j^. Then the graphs A^_j^, i G {a,b,c} are glued up as 
described above. Thereafter, three decorating bonds are drawn to connect the external nodes. 
This procedure is repeated ad inhnitum. 

2.2 Definitions 

In this subsection we begin performing the mathematical construction of the model outlined 
above. In order to fix the terminology, we recall relevant mathematical notions. A simple graph 
G is a pair of sets (V, E), where V is the set of nodes, whereas E is a subset of the Cartesian 
product V X V. It is symmetric and irreflexive, i.e., {j,i) G E whenever {i,j) G E, and {i,i) ^ E 
for every i,j G V. We say that i and j are connected by a bond if {i,j) G E. In this case, we 
write i ^ j and say that i and j are adjacent or that they are neighbors. Hence, the elements 
of E themselves can be called bonds. The graph is said to be complete, if each two nodes are 
adjacent. For a given i, by n(f) we denote the degree of f - the number of its neighbors. If 
V, and hence E, are finite, the graph is said to be finite. Otherwise, the graph is infinite. An 
infinite graph is called locally finite, if n{i) is finite for every node. 

Given G = (V, E) and G' = (V', E'), let </> : V —V' be such that (/>(i) ~ (/>(j) whenever i ~ j. 
Such a map cj is called a morphism. A bijective morphism is called an isomorphism. If cj is 
an isomorphism, then its inverse is also an isomorphism, and then the graphs G and G' 
are said to be isomorphic. Such graphs have identical structures and thus can be identified. In 
this case, we also say that G' is a copy of G. One observes that this refers to both finite and 
infinite graphs. An isomorphism iji : V —)• V, i.e. which maps the graph onto itself, is called an 
automorphism. The graph G' = (V', E') such that V' C V and E' C E is said to be a subgraph of 
G = (V, E). In this case, we write G' C G. Suppose that a subgraph G' C G has a copy, say G", 


2 












that is, there exists an isomorphism (f) : G" ^ G'. Then cf), considered as a map 0 : G" —)• G, is 
called an embedding of G'' into G, whereas G' is called the image of G" under this embedding. 
Fig. [1] presents the so called unlabeled graphs, which are studied in this work. After labeling, 
i.e., attaching a label to each of the nodes, such a pattern turns into a graph. Another labeling 
may or may not give the same graph up to an automorphism. This depends on whether or not 
there exists the corresponding automorphism. 



Figure 3: Labelled graphs of level one based on motifs Mi, M 2 , M 3 , M 4 , M 5 . 


For instance, any labeling of the triangle Mi yields the same graph as in any case each of the 
nodes has the same neighbors. So the triangle has six automorphisms. For the pattern M 2 , the 
corresponding graph shown in Fig. [3] with the interchanged labels a and b is the same. However, 
the graph with the interchanged c and d is not the same anymore. Of course, this new graph is 
isomorphic to the initial one. This is because there is only one nontrivial automorphism of M 2 '. 
the one which interchanges a and b, and preserves c and d. 

Let G' C G and G" C G and there exists an isomorphism (/> : G' —)• G". Then we can consider (j) 
as an equivalence G' ~ G". The equivalence class of G' is defined as the set [G'] = {G" : G' ~ G"}. 
It is called motif. If the number of appearances of motif [G'] in a network G is higher than the 
number of its appearances in the Erdos-Renyi random graph, then [G'] is called network motif. 

Now we present the notion of a random graph, which we use in this work. The random graph 
model is dehned to be a pair consisting of an underlying graph G = (V, E) and a probability 
space (E,T,P). If G is finite, as £ one can take the set of all subsets of E. In the sequel, we 
deal with such random graph models only. Thus, for E' E T, we say that E' has been picked 
at random with probability E(E'). In many models, the bonds are being picked independently 
with probability which may depend on the bond. In this case, one deals with a random graph 
model with independent bonds. For such graphs, 

P(E') = n p(e), (2.1) 

eCE' 

where p{e) is the probability of picking bond e. The set of graphs 

{G' = (V,E')}E'ef 

is called the graph ensemble - each G' is being picked at random from this ensemble. Now 
suppose that we have two random graph models with independent bonds. We have to specify 
the definition of isomorphism for such graphs. Let Gi = (Vi,Ei) and G 2 = (V 2 ,E 2 ) be their 
underlying graphs and pi,P 2 be their corresponding probability (|2.1I) . Then the map (/> : Vi —V 2 
is said to be the isomorphism of the random graphs if there exists isomorphism / : Gi —)■ G 2 (in 
the meaning shown previously for non-random graphs) such, that for every {i,j) E E we have 

Piiihj)) =P2 {{f{i),f{j))). 


2.3 The construction 

As was mentioned above, each of our graphs is constructed in an algorithmic way from the 
corresponding motif presented in Fig. [TJ As they are random graphs with independent bonds, 
we have to construct the corresponding underlying graphs and to define the probability of picking 
the bonds, cf. (|2.1I) . In all our models, the bonds will be of two kinds, which we call basic bonds 
and decorations. Basic bonds are non-random, i.e., picked with probability one. Decorating 
bonds appear with probability p E [0,1], which is a parameter of the model. Now we present the 
formal construction of the underlying graphs. Let q be the number of nodes in the corresponding 
motif, i.e., q = 3 for Mi and g = 4 for the remaining motifs. At step k = 1, we just label the 
nodes of the corresponding motif hy i = 1,..., q and obtain the initial graph Ai = {Vi,Ei). All 
its bonds are set to be basic. Suppose now that we have q + 1 copies of Ai obtained by the 
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isomorphisms j = 0,1,... , g. Thus, in j-th copy the nodes are i = 1,... ,q. The graph 

A 2 is obtained from these copies under the following conditions 

* = 02 (j) = '/' 2 (*)> * = ( 2 - 2 ) 

Thus, the images of Vi under 4>2 and 02 with i ^ j intersect only at one node where (|2.2p holds. 
The maps 02 ! J = 0i 1; • • • i 9 embed Ai into A 2 . The nodes 02 (Oi * = li ■ ■ ■ i O'! are called the 
external nodes of A 2 . All other nodes are called internal. Thus, A 2 has q external and q{q — l)/2 
internal nodes. At this stage, we label them by i = 1 ,... ,q{q + l )/2 in such a way that the 
external nodes have the same labels as in Ai, that is, 02 (i) = i, i = 1,... q. By construction, the 
bonds obtained as images under the map 02 are decorations: they are of the form ( 02 (Oi 
where i and j are adjacent in Ai. From the first condition in (j2.2l) we see that the decorating 
bonds connect the external nodes of A 2 . The remaining bonds of A 2 are set to be basic. Now 
we construct A*, for A: > 3 from one copy of Ai and q copies of Afc_i. Let 0° be the map which 
produces the copy of Ai and 0 ^, j = 1 ,... be the maps which produce the copies of Afc_i. 
We then impose the conditions 

* = 1, •••!?; = </'i(*)! * = l,- • • !^! * / j (2-3) 

and obtain A^. Thus, 0° embeds Ai —>■ A^, and 0^ : Afc_i —>■ A^., i = 1, 2,... , g. As above, 
the nodes 0^(i) are set to be external, and the remaining nodes are internal. The images of V 2 
under 0^ and 0^ with i ^ j intersect only at one node where (j2.3l) holds. Again we label the 
nodes of A^, is such a way that 0^(i) = i, i = 1,... ,q. Now let us establish which bonds of 
Afc_i are decorating and which are basic. As above, the bonds connecting the external nodes 
are decorating. The images of decorating bonds of Afc_i are decorating bonds in A^.; the same is 
true also for the basic bonds - the basic bonds of A^ are exactly the images of the basic bonds of 
Afc_i. As above, by I 4 and we denote the sets of nodes and bonds of A^, respectively. Thus, 
for A > 2 we have Ek = E'l^U E'^, where E'^ (respectively, E'^) consists of basic (respectively, 
decorating) bonds. All A^, A: E N, are considered as subgraphs of an infinite graph A^o, the 
structure and properties of which are not important for the study presented in this article. 

Note that the construction principle used above essentially differs from that used in PEI 
El Eg. Namely, in our case to obtain A^ one replaces each node of the basic pattern by a copy 
of the graph Afc_i. In the hierarchical lattices, one replaces a bond. As we shall see in the sequel, 
this leads to essentially different properties of the resulting graphs. Below in Fig. m we illustrate 
the construction described above for the case where the basic pattern is the motif Mi. In this 
case, the bare graph (which occurs for p = 0) is the approximating graph for the Sierpihski 
triangle. The elements of E 2 (middle graph) and of E'^ (right-hand graph) are depicted as solid 
lines, whereas the elements of E 2 and of E'^ appear as dotted lines. We omit some dotted lines to 
indicate that they are random and hence may be absent in a given realization of the graph. Note 
that A 3 can be viewed as the triangle composed from three copies of A 2 . In Fig. [Sj we present 
the construction of the bare graph A 3 corresponding to M 2 . In contrast to the former case, it is 
not a planar graph. In Fig. 4, we construct the bare graph A 2 for motif M 3 . One observes that 
in this picture the node c of the lower left-hand quadrat (i.e. quadrat o) is glued up with node a 
of the upper right-hand quadrat. It is interesting that the corresponding fractal can be obtained 
by the following procedure, resembling the one which yields the Sierpihski triangle. One takes 
the full quadrat and cuts it out into four equal quadrates, not cutting the external lines. Then 
one glues up the vertices of the smaller quadrates as depicted and proceeds with cutting out the 
smaller quadrates. The fractal which one obtains from M 5 is a three dimensional version of the 
Sierpihski triangle. One takes the full tetrahedra and cuts out its inner one fourth in such a way 
that the remaining four tetrahedrae are glued up according to the rule: vertex b of tetrahedra a 
is glued up with vertex a of tetrahedra b, etc. 

2.4 Degree distribution 

Now we turn to the structural properties of the graphs constructed above. Here and subsequently, 
q and r stand for the number of nodes and bonds in the corresponding motif, respectively. By 
the construction described above, the number of nodes in A^ is |I4| = ( 7 |I 4 _i| — q{q — I)/2 and 
|IA| = q. Likewise, the expected number of bonds is I Fill = r and {\Ek\) = q{\Ek-i\)+ rp. Hence 

\yk\ = ^^, m\)=rq^-^+rp^^-^, A E N. (2.4) 
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Figure 4: Construction of the graph A 3 based on Mi 




Figure 5: Construction of the bare graph A 3 based on M 2 


As was mentioned above, the degree distribution is an important characteristic of the graph. 
In contrast to the Erdos-Renyi type graphs, the distribution of the random variable n{i) in our 
graphs depends on the type of i. Thus, the simplest way to describe this distribution is to 
average n{i) over the nodes of a given A^, that is, to consider Uk = |I 4 |~^ ^(*)- Let (uk) 

be the expected value of Then 


(nk) = 2 


mi) 

\Vk\ 


4r 

9(9-1) 



9 - 1 + 2p \ 


(2.5) 


However, this result gives only partial information about the node degree distribution. To 
get more let us analyze the structure of the node sets I 4 , k = 1,2,.... For a given A^ and 
I = let be the set of nodes i £ Vk which have the same degree distribution, 

independent of k for I < k — 1. For the graphs based on regular motifs Mi, M 3 and M 5 , 
consists of the nodes which are external for some A/ and, at the same time, are internal for any 
A;_|_i. Here we mean those A;’s which are subgraphs for A^. As an example, let us consider the 
graph A 2 based on Mi, see the middle graph in Fig. HI The nodes a, b, and c constitute V 2 ^ 




Figure 6 : Construction of the bare graph A 2 based on M 3 
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whereas the remaining nodes constitute V 2 ^\ For M 2 and M 4 , this partition is more complicated 
and will be discussed below. First we analyze based on complete motifs Mi and M 5 . The 
elements of are exactly the nodes at which the subgraphs j = 1 ,..., g are glued up 

to form A/^, whereas the elements of ^ are exactly the external nodes of A^. Then = q 

and = q{q — l)/2. For I < k — 1, we have = g|14^i\l, which can be solved to yield 

\V^'^\ = ^q^-\q-l), l = |yf^|=g. ( 2 . 6 ) 

The degrees of i G are non-random as these nodes receive no decorating bonds. For such 
i, n{i) = Xj where n^^\j) is the degree of the corresponding node in the basic pattern, 

and the sum is taken over all such patterns which are glued up. By the symmetry of Mi and 
M 5 , we have that n(i) = 4 for Mi and n{i) = 6 for M 5 . For i G / = 2, 3,... , fe — 1, we have 
n{i) = n{i) + ^{i), where h{i) is non-random and has to be calculated as just described. The 
summand n(z) is the number of decorating bonds attached to i. For / = 1,... , A: — 1 and i G 
we have h{i) = 2{q — 1 ) and v{i) takes values n = 0 , 1 , 2 ,... , 2 (g — !)(/ — 1 ), with probability 

Prob = n) = (2.7) 


For i G v{i) takes values 0,1,... , (g — l){k — 1). As is usual in the theory of real world 

networks, which are in fact non-random, the randomness manifests itself as the random choice 
of a node. If we apply this principle here, then (12.711 can be considered as the conditional 
probability distribution, conditioned at the event that the node i is been picked from the set 
The probability of this event is taken to be proportional to the number of elements, that 
is, 


Prob (^i G 
Prob (i G 

Now we take the expectation of n{i) 

(nk) 


\Vk\ 


g- 1 
1 + q^~^ 



I < k — 1, 


2 

gfc-l _|_ I • 

with respect to this distribution and obtair0 
_ q’^-^{2q-2 + 2p) - 2p 

~ g^-i -b 1 ’ 


which agrees with (12.51) . In the same way we find the second moment 


( 2 . 8 ) 


(2.9) 


(n^) = 4(g - 1)^ (8g - 6)p (4g -b 2)p^. 


In order to hgure out the limit k —>■ -|-oo of the distribution given by (1231) and (12.8p we calculate 
its characteristic function, cf. 


k-l2{q-l){l-l) 


^k{t) = Y. Y. 


1=1 u =0 


J2(q-l)+u)it i^ 1) • g 
1 -b q^-’^ 


^2(g -!)(/- _ p-^2{q-l)il-l)-u 


u=0 

{q-l){k-iy 


g^-1 -b 1 


Then the limiting characteristic function is 

(g - l)e2(‘'-i)i* 


q-l){k-l)-L/ 


p{t) = 


q — (e**p -b 1 — 


^Detailed calculations of this and similar quantities are given in Appendix below. 


( 2 . 10 ) 
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As was mentioned above for the graph based on M 3 , the same node partition can be applied also 
for the graphs based on Mi and M 5 . Motif M 3 is regular hence (|2.6p and (I2.8p still holds. Here for 
I = 1,2, —1 and i G we have n(i) = 4+z^(i) and u{i) takes values = 0,1, 2,... , 4(Z — 1) 

with probability 

Prob(i/(Z) = u)= 

The external node i € has random degree ^{i) = 0,1,..., 2{k — 1). Taking the expectation 
of n{i) with respect to this distribution one obtains, see Appendix, 


(nk) = 4 + 


p- 


3 + 2p \ 
4^-1 + ly ’ 


( 2 . 11 ) 


and the k —+00 limit of the second moment 

(n^) = 16 + 12p + —p 


68 2 
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The characteristic function is 
^k{t) = 


3e 


Ait 


+ 


1 + 4i-fc 
2 g 2 it 


1 _ 4I (giip _|_ 1 

4 — + 1 — pY 


A{k-1) 


+ 


which in the limit k 


4fc-i 1 

00 takes the form 
<p{t) = 


(e 'V + 1 - p) 


2{k-l) 


= V^, 


3e 


4it 


( 2 . 12 ) 


4 — (e^*p + 1 — p)^ 

It remains to analyze the graphs based on M 2 and M 4 . Label the nodes of M 2 as shown in the 
Fig. O For the external nodes of the corresponding graph, we have 


n(a) = n( 6 ) = 2 + zz(a), zz = 0 , 1 ,..., 2{k — 1 ), 

n(c) = 3 + zz(u), zz = 0,1,..., 3(Z — 1), 

n{d) = 1 + i^{d), z^ = 0,1,..., Zc — 1. 

For each I = l,2,...,/c — 1 set the consists of three subsets of the same cardinality with the 
following degrees 


This yields, see Appendix, 


n(i) = 3 + zz(i), zz = 0,1,... , 3(Z — 1), 

n(i) = 4 + zz(i), zz = 0,1,... ,4(/- 1), 

n(i) = 5 + zz(i), zz = 0,1,... , 5(Z — 1). 


and the second moment 

For the characteristic function, we have 


, 50 112 214 2 


Pkit) = 


1 


1 + 41 -*^ L 


o 3 it 


1 _ 4 I _ p^3(fc 1 ) 

4 — {pe^^ + 1 — p)^ 


+ 


+ 


+ 

+ 


1 - 41 -^(pc'* + 1 - p)4(*^-l) 

4 — (pe^* + 1 — p)^ 

1 - 41 -*^ (pA* + 1 - p)5(*=-i) ] 
4 — (pe^* + 1 — p)^ 

2 e2'*(pA' + l-p)2 (^-^) + 


+ 


4^ + 4 

g3it(pgii ^ ^ git(^git ^ 


(2.13) 
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which in the limit /c —>■ oo yields 


oSit 


Ait 


/ \ ^ ^ 

^-Q H -^- A 

4 — (e'*p + 1 — p) 4 — + 1 — p) 

g5it 

+ 


+ 


4 — (e^*p + 1 — p) 

For the graphs based on M 4 , we obtain 


5 • 


m=5+|(p-jI4^) 


, 2 ^ 76 167 335 n 

in ) =- 1 - p-\ -» , 

' ' 3 9 27 


(2.14) 


p{t) 


Mt 


+ 


2 e 


bit 


2 (4 — + 1 — p) ) 4 — (e'*p + 1 — p) 


+ 


+ 


g6ii 

2(4 — + 1 — p)® ) 


(2.15) 


For all our graphs, the limiting characteristics functions can be continued to functions analytic 
in some complex neighborhood of the point t = 0. This means that the limiting node degree 
distribution has all moments and hence cannot be of scale-free typ^. Another observation here 
is that the characteristic function of the Poisson distribution 


7?Poisson(t) = exp [c (e‘* - l)] , 

can be continued to a function analytic on the whole complex plane. Therefore, the degree 
distributions in our graphs with p > 0 are intermediate as compared to the Poisson and scale- 
free distributions. For p = 0, our functions (I2.10j> . (I2.12j> . (I2.14j> . (I2.15p are also entire. 


2.5 Amenability 


The property of our graphs which we address now is amenability. Let G = (14, E) be a countable 
graph with node set V and bond set E. For a finite Z\ C F, by dA we denote the set of nodes 
which are not in A but have neighbors in A. By | A| and \dA\ we denote the number of elements 
in these sets. The graph G is said to be amenable if there exists a sequence of finite node sets 
{Ak}ken, such that 


lim 

fc^ + OO 


\dAk 

\^k\ 


= 0 . 


(2.16) 


If such a limit is positive for any sequence the graph is called nonamenabile. Some¬ 

times, sequences for which ()2.16l) holds are called Van Howe sequences. Cayley trees, except 
for Z, are nonamenable. Let us consider the underlying graphs of our random graphs. Due to 
their hierarchical structure, it is convenient to check (|2.16H for the sequence of node sets of A^, 
that is for { 14 }fcgN- by the construction of A^., the inner boundary of each Vk is the set of all 
its external nodes, the number of which is equal to the number of nodes in the corresponding 
motif, i.e. it is q. By construction, g — 1 of them become inner nodes of Afc_|_i, and receive new 
k{q — 1) neighbors (and none in the next steps). The remaining one becomes an external node 
of Afc_|_i. We can choose {FfcjfcgN in the way that this external node becomes an internal one 
in the next step. And then it receives new {q — 1)(A; -|- 2) neighbors outside A^ and none in the 
next steps. Then for all the graphs we obtain 


I^Afcl _ k{q - 1)^ (g - l)(fc 2) fc^oo 

|Afc| U^’^ + q ) ^ 


which means that all our random graphs are amenable with probability one. 

^For scale-free graphs, the node degree distribution is P{k) = Ck fc > 1, 7 > 1; hence, 'YAk=ik^P{k) 
diverges for all m > 7 — 1 . 


























2.6 Clustering 

For a given node i ^ V of degree n(i), let N[i) be the number of bonds linking its neighbors 
with each other, which is the number of triangles with vertex i. Clearly, N{i) < n{i)[n{i) — l]/2 
and the maximum value of this parameter is attained for complete graphs where each node is a 
neighbor to all other ones. Thus, the quantity 


Q(t) := 


2N{i) 


n{i)[n{i) — 1 ] 

characterizes clustering at node i. Then the clustering of our graphs we define as 


Q 


lim 

k—^+oo I I 


ieVk 


Note that for many graphs, e.g., for trees or bipartite graphs, one has Q{i) = 0 for any node i, 
see also [MlE?]. For random graphs, the degree n(f), as well as the parameter N(i), are random. 
The calculation of Q in this case is much more involved. We will address it in a forthcoming 
paper. Here we only compare the values of Q obtained for the bare and fully decorated versions 
of our graphs, i.e., for p = 0 and p = 1. 

For the bare graph based on Mi, we have n{i) = 4 for internal node f G 14 n{i) = 2 
for external node i G Vk- Besides 


N{i) 


3 

< 2 
1 

\ 


iev. 


( 1 ) 

■ C ’ 
r{l) 


i G 1 = 2,3,. 


i G H 


(k) 


k — 1, 


which follows directly from the construction of the graphs. By 


(12.411 and (12. 6 p one gets 


1 


ieVk 


1 


Z] Z] = 

l = l 


I l^.l-I^Zl-3 3 1 2 

2|44| 3|44| iFfcl 3 6|44| iVfcl' 


Hence, 


Q 


lim 

k^-\-oo I V}^\ 


Z *3(0 

i&Vk 


1 = 0 ,4444... . 


For the fully decorated graph based on Mi, internal node i G I = 1,2,... , k — 1, has degree 

n{i) = 4/ and N{i) = 4/ whereas for external node i G we have n(i) = 2k and N{i) = 2k —1. 
Then 


1 


Z *3(0 

i£Vk 


1 

m 


Z Z = 


1 = 1 


1 


fv V 2 ( 2 fc-l) ^ 

41(4/-1) ' ^ '2k{2k-l) 


1 


2 • 3' 


fc-i 

E 

1=1 


3"^ 
41 - 1 



4-3*^ ^ 3-^ 6 

3^ + 3 ^ 41 - Z A:(3^ + 3)' 


Hence 


(5 = 2-3 arctan 3 


3-1/" In 


3"/" + 1 
31/4 _ 1 


0,525897. 


For the bare graph based on M 5 , one obtains for internal node i G V^^\l = 1,2,... ,k — 1: 
n{i) = 6 , N{i) = 8 for 1 = 1, and n{i) = 6, N{i) = 6 for 1 > 2. For the fully decorated graph 
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based on M 5 we have for all internal nodes n(i) = 61, and N{i) = 8 for I = 1, and N(i) = 121 — 3 
for I > 2. Hence, for the bare graph, we get 

/2 2\Vu^\ 12 \ 

and for the fully decorated graph 

Q « 0.554145, 

which surprisingly is quite close to the clustering in the bare version (|2.17l) . 

2.7 Small-world property 

There exists one more property of real-world networks which Erdos-Renyi type graphs do not 
share, see e.g. [T8l[T9]. It is the so called small-world property. To formulate it one needs the 
following notion. A path in the graph is a sequence of nodes such that every two consecutive 
elements are neighbors to each other. The length of the path is the number of such consecutive 
pairs, which is equal to the number of bonds one passes on the way from the origin to the 
terminus. If every two nodes can be connected by a path, the graph is said to be connected. 
For a given two nodes, i and j, the length of the shortest path connecting them is said to be 
the distance p{i,j) between these nodes. Informally speaking, a graph G = (y,E) has the 
small-world property (is a small-world graph) if every two nodes i,j&V are ‘not too far’ from 
each other. More precisely this property is formulated as follows. An inhnite graph G has a 
small-world property if there exists a sequence of its connected hnite subgraphs with 

the following property. Let diam(G'fc) = maxjjgVfe p{hj) be the diameter of Gk, A: G N, and (ufe) 
be the average value of the node degree in Gk, that is, (n^) = 2|Efc|/|I4|. Then the sequence 
{GfcjfcgN) and hence the graph G, are said to have the small-world property if there exists a 
positive constant G such that for all A; G N 


diam(Gfc) < Clog^^^) |I 4 |- 


In such graphs, the distances between the nodes scale at most logarithmically with the size of 
the graph. Let us consider this characteristic of our graphs without decorations, i.e., for p = 0. 
For a chosen motif, the diameter or A^ is the maximum distance between two external nodes 

diam(Afe) = maxp(i,j), i,j G / j. 

For the complete motifs Mi and M 5 there is diam(Ai) = 1, and for the other motifs diam(Ai) = 
2. By the construction of A^, k = 2,3,..., it is easily seen that the distances between two 
external nodes increases two times ay each step. Hence 

diam(Afc) = for Mi and M 5 , 

diam(Afc) = 2^ for M 2 , M 3 and M 4 , 

that means that the diameters scale exponentially with the size of the graph. 

For p = 1, the distance between two chosen external nodes in A^, k = 1,2,..., is 1. Hence, 
we have to analyze the distances between other pairs of nodes. Here we present the results for 
graphs based on motif Mi only. The distance between an internal and an external node in A 2 
does not exceed 2. The distance between two such nodes in A 3 is not greater than 3, and k in 
Afc. Therefore, to estimate the distance between two internal nodes in Ak, one has to find the 
greatest / < A; — 1 such that these nodes belong to different A/. Then add the distances between 
these nodes and the common external node. Hence 

diam(Afc) < 2{k — 1). 

Thus, neither of our bare graphs has the small-world property. At the same time, this property 
holds for all fully decorated graphs. 
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3 Phase Transitions in the Ising Model 


There exists a deep connection between the properties of Gibbs random fields of the Ising model 
and the structural properties of the underlying graphs, see [6]. In the physical terminology, 
each (pure) Gibbs random field corresponds to a state of thermal equilibrium of the model, see 
[20] for more details. Accordingly, the existence of multiple Gibbs random fields corresponds to 
the existence of multiple equilibrium states and hence to phase transitions. For noninteracting 
spins, the Gibbs random field is unique. However, if the interaction is strong enough and if it 
is effectively propagated by the underlying graph (due to high ’connectivity’), the Gibbs fields 
can be multiple. 

The Ising model on an infinite graph G = (V, E) is defined by assigning spin variables 
CTj = ±1, i G F. Two spins, crj and Gj, interact whenever i ^ j. The space of spin configurations 
is then S := {—1,1}^. It is equipped with the discrete topology and the corresponding Borel 
cr-field. A Gibbs random field is defined as a probability measure on S which satisfies a certain 
condition formulated by means the so called Gibbs specihcation, see [20]. The specihcation in 
turn is constructed by means of conditional model Hamiltonians, defined as follows. For a finite 
A C F, ^ G S, and a fixed value of the inverse temperature /3 > 0, the conditional Hamiltonian 
in A is given by the following expression 

- /3'Hz\(crz\|^) = ^ + ^ ^ (3.1) 

where = {ai : i G A}, h is an external held, and Jtj G M is the spin-spin interaction intensity. 
Note that the latter parameters include /3. For hierarchical graphs constructed in an algorithmic 
way, the inhnite graph (F, E) is obtained as limiting object, dehned by means of a system of 
embeddings which map each hnite fragment into the graph, see mm- In this article, however, 
we do not follow this way and consider the annealed case, see m, where one deals exclusively 
with states on such hnite fragments. Here we mean the randomness related to the decorating 
bonds. 

In view of the hierarchical structure of our graphs, we take A = A^ in (13.111 to be the set 
of the inner nodes of a given A^ corresponding to motif Mi. Then j in the last term in (|3.1I) 
runs through the set of external vertices of A^. We also restrict our consideration to the case of 
h = 0. For A: = 1 we have no internal nodes and no randomness either. Thus the corresponding 
Hamiltonian is 

— iSTii = K{ab -|- oc -|- be), 

here we use the shorthand like a = Ga, and K stands for the interaction intensity corresponding 
to nonrandom bonds. Recall that, for A; > 1, by E'f, (resp. E'j^) we denote the set of solid (resp. 
decorating) bonds of A^. To take the latter randomness into account we introduce independent 
random variables u G {0, such that Prob(a;ij = 1) = p. Then we set Jij = = K for 

{i,j) G E'^, and Jij = = Lujij for (i,j) G E'j^. In general, we assume that K ^ L as the 

random and nonrandom bonds play different roles in our constructions. The nonrandom bonds 
form a skeleton of the graphs, whereas the random ones increase its connectivity. Moreover, by 
setting K = 0 we can pass to the model defined on a purely random graph. 

Then, for A: > 2, (|3.1I) takes the form 

-P'Hk{crAt\a,b,c) := -/3^zifc(o'zifc|a, 6, c) = L {uatab + ujacac + ujbcbc) (3.2) 

-1-0 ^ Jai^i + b ^ J^Gi -I- C ^ J^CTj -I- ^ JfjGiGj, 

ieNl 

where NJ! = {i G F^ : i ~ u} is the set of the neighbors of v in A^ and E^^ C E^ is the set of 
bonds connecting the inner nodes of A^ to each other. 

The Hamiltonian in (13.2p can be rewritten in a recursion way, see Fig. [T] as follows 

-l37ik{GAk\a,b,c) = L{u}abab + ujacac + u}bcbc) - l37ik-i{GAi ,\a,'y,l3) 

(3.3) 

- P'Hk-i{G^b^_J-f,b,a) - l3'Hk-i{GAi_J/3,a,c), 
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Figure 7: Graph A^, 


where stands for the set of inner nodes of x G {a,b,c}. Then the annealed Gibbs 

distribution in is defined as follows 

T^ki^Aja,b,c) = -—^-—exp{-l3'Hk{(JAja,b,c)), k>2, 

Z/kyCl-! 0, C) 

and TT^ = 1 as Ai = 0. Here the partition function has the form 


Zk{a,b,c) = exp {-fj'HkjcrAk\a-,b,c))\ 


Zi{a, b, c) = exp {K{ab + ac + be )), 
and (•) denotes the expectation in uj. 

For A: > 2 let / : {—1,1}^'= ^ M be a local observable, which is a function dependent on 
with some m < k such that Am C A^. Set 


Fk{f\a,b,c) = (^fi(^Ak)'^ti(^Aja,b,c)^, k > m. 


(3.4) 


The sequence {Fk{f\a,b,c)}k>m is bounded and thus has accumulation points. Our aim is to 
study their dependence on the values of the boundary spins a, b, c. 

In view of the independence of the bond variables w, we have 

Zk{a,b,c) = ^ exp (La;a6a5) ^ ^ exp (Lwacflc) ^ ^ exp (Lwfecbc) ^ (3-5) 

X ^ Zfc_i(a,7,/3)Zfc_i(7,6,a)Zfc_i(/3,a,c). 

Assume now that the observable / depends on the spins indexed by Am C Then by (|3.3I) 

and (|3.4I) we have 

^ ^ exp (LuJabab) ^ ^ exp (Lujacac) ^ ^ exp (Lubcbc) ^ x 


Ffc(/|a, 6 ,c) = 


E ( E f{(TAi_Jexp (^-/3'Hfe_i((T^«_Ja,7, 

o,/3,7 

exp (^-/377fc_i(<T^6_j7,6,a)^ ^ x 


^ exp ( -13'Hk-i 1 ^, a, c 




k > m, 


that can be rewritten as follows 
Fk{f\a,b,c) = 

_ Ea,/ 3,7 7, l3)Zk-i{a, 7 , j3)Zk-i ( 7 , b, a)Zk-i (/3, a, c) 

T.a,f}a Zk-i{a, 7 , f3)Zk-i ( 7 , b, a)Zk-i {13, a, c) 


(3.6) 
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From now on we assume that the locality of / is such that it m = 1, see dm), and that the 
corresponding Ai is a subset of for all k > 2. Then, in addition to (|3.4I1 . we set 

Fi {f\a, b, c) = f{a, b, c) = f{a, c, b) > 0, (3.7) 


where we assume also that / is positive and symmetric with respect to 6 •H' c. Next we introduce 
the following variables 


M '■= 

Bk '■= 

yk ■= 
Yk ■■= 
YS := 


Ffc(/|l,l,l) 

Ffc(/|l,-l,-l) =Tfc(/|-1,1,1), 


and 


= (pe^ + 1 — p) ^ , 

R~ = (pe^ + 1 — p) {pe~^ + 1 — p)^ 
Then from (j3.5p and (j3.6p we get 

Afc+i = R+{Al + iAkBl + ^Bl), 
Bfc+i = R-[Al + ^AkBl + Wl), 


and 




fc+i 


fc+l — 


fc+l 


YjAkjAl + Bl) + 2Y^Bl{Ak + Bk) + 2Y,^Bl 
Al + MkBl + 4Bl 

Y^AkBkjAk + Bk) + 2Y^Bl{Ak + Bk) + Y^B^Ak + Bk) 
BkAl + AAkBl + Wl 

2Y^AkBl + 2Y^Bl{Ak + Bk) + Y^Bk{Al + Bl) 


BkAl + 4AkBl + Wl 


with the initial conditions 

and 

Then for 
by (I3.9P we get 
where 

and 


Al = e' 


3K 


Bl = e 


-K 


Yl=f (1,1,1), Yl = fil,l,-l), Yl = fil,-l,-l). 

Xk = ^k/Bk, xi = e'^^>0, 

Xk+l — t4^(xk)j 


(pix) = 


+ 3x + 4 — X + 4 

x^ + 4x + 3 


t = 


RA 


pe^ + 1 — p 
pe~^ + 1 — p 


X + 3 
2 


In these notations, ()3.10p can be rewritten in the following form 

Ifc+i = T(xk)Yk, /c E N, 

where Yk is the column vector transposed to (Y^,Yl,Yl) and 


T(x) = 


x(x^ + l) 

2(x+l) 

2 

x^+3x+4 

x^+3x+4 

a:^+3a:+4 

x{x+l) 

2(x+l) 

x+1 

qF+4x+3 

x^+4x+3 

aT+4x+3 

2x 

2(x+l) 

x'^+l 

x^+4x+3 

a:^+4a:+3 

x^+4S+3 


(3.8) 


(3.9) 

(3.10) 


(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 


13 





















Observe that, for each x > 0, T{x) is a stochastic matrix, which means that each of its rows 
consists of nonnegative elements and sums up to one. Then, for each k G'N, the matrix 


Sk = T{xk)T{xk-i) ■ ■ ■ T{x 2 )T{xi) 
is also stochastic, and the solution of the recursion in (|3.15ll is 

Yk+i = SkY^. 


(3.17) 


(3.18) 


Products of stochastic matrices as in (I3.17p appear in the theory of inhomogeneous Markov 
chains, see e.g. |22l [23l I24j . They also are being used in communication networks, control 
theory, parallel computing, and decision making, see [251 ESI Ea EH] and the references therein. 

As mentioned above, our aim is to study the limits of the sequences {T)!}, i = 1,2,3, 
defined in (13. 8p - (|3.1ip . and hence described by (|3.15p . (|3.18p . If, for an arbitrary Yi, the 
limit Too = hmfc_>.oo is a vector with all components equal to each other, then the limiting 
average (1331) is independent of the boundary spins, which corresponds to the uniqueness of 
the limiting Gibbs state, and hence of the state of thermal equilibrium of the model. In the 
terminology of Markov chains, this is related to the ergodicity of the sequence {T(xfc)}. By 
definition, see e.g. |28l Definition 1, page 1479], such a sequence is ergodic if the product 
sequence {Sk} as in (|3.17l) converges to a stochastic matrix with identical rows. In [26], such a 
sequence is called consensus. In this case, the sequence {Yk} converges to a vector with identical 
entries. Likewise, the existence of subsequences {Y^^/eN convergent to vectors with nonequal 
components corresponds to the multiplicity of such states, and hence of a phase transition. Since 
we do not introduce the Gibbs states of our model explicitly, we use the following 


Definition 3.1. For fixed K, L i p, the Ising model on our graph is said to be in an unordered 
state if for each observable / satisfying dSZD, there exists linik^+oo Fkif\o-,b,c) independent of 
b and c, and hence of a. Otherwise, the model is said to be in an ordered state. 


As might be seen from (|3.15p . the limiting properties of the sequence are crucial for 

the corresponding properties of {Yfcjfcgp^. 


Lemma 3.2. Assume that the sequence {xfcjfceN defined in h3.12\) - {3.14^ converges to a certain 
X* > 0. Then the sequence {T{xk)}k&N defined in 113.16\} is ergodic. 


Proof. By definition, each row = (th, ..., Tin) of a stochastic matrix T = {rifinxn is a proba¬ 
bility distribution. For two such rows, we define 

1 ” 

IIA — TjW = - \Til — Tjl\. 

1=1 

Then the Dobrushin ergodicity coefficient of T is 


D{T) := max Ht, — r,-||. 

It can also be written in the form; 

n 

D{T) = 1 - miny^ mm{Tii]Tji}, 


which yields that D{T) < 1 whenever all Tij are strictly positive, see [27] for more detail on 
this issue. On the other hand, for two stochastic matrices T and Q, it is known that, see |23l 
Ghapter 3], 

D{TQ) < D{T){Q). (3.19) 

Since the matrix elements of T(x) in (|3.16p are continuous in x > 0, one has T{xk) —>• T(x*), 
component-wise, as k ^ -|-oo. Therefore, each element of the latter matrix is strictly positive, 
which yields D[T[xfij) =: 5 < 1. By the mentioned continuity we also have that, for a given 
e > 0 such that <5 -|- e < 1, there exists k^ such that D{T{xk)) < 5 -|- e < 1 for all A: > This 
yields by (I3.17p and (I3.19p that D{Sk) —)• 0. As a sequence of stochastic matrices, {Sk}keN 
contains convergent subsequences, each of which converges to a stochastic matrix S with strictly 
positive elements, for which D{S) = 0. By the Perron-Frobenius theorem there exists only one 
such limit. □ 
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For a given k gN, let d{Yk) be the diameter of Y^, that is, 

d{Yk) ■■= max (Y^-Yj!). 
i,j=l,2X 

It is known that, see m, for each Yi, 

d{SuYi) < D{Su)d{Yi), 

which yields d(Yfc) —0 as A; —oo if {T{xk)}keN is ergodic. On the other hand, for each e > 0, 
one finds Yi with positive entries such that 


d{SkYi) > [D{Sk)-e]d{Y^). 


(3.20) 


As possible limits of {xfcjfcgjs}, there can appear the solutions of the following equation: 


X = t(j){x), X > 0. 


(3.21) 


First, let us consider the case of L < 0 where the interaction along the random bonds is antifer¬ 
romagnetic. Then t < 1, see (|3.14l) . and the only solution of (|3.21l) is 


-3-t + V9 -h 22t - 15t2 

^ ’ 


which is clearly positive. By (|3.13|1 we have 


(3.22) 


Therefore, the solution (j3.22D is stable and Xk x* for all xi > 0. Thus, by Lemma 13.21 the 
model is in an unordered state for such L and all AT G M and p G [0,1]. This possibly holds 
due to the frustration caused by the motif the graph is based on. For L = 0, i.e., for the graph 
without decorations, we have t = 1 and the only solution of (|3.21l) is x* = 1. In this case, the 
model is in an unordered state for all K gM. 

Let us turn now to the case of L > 0, in which t > 1. Now (I3.2ip has two solutions 


^(1) _3 + t- V9 22t - 15t2 ^(2) _3 + t + V9 -k 22t - 15t^ 

~ 2(r^i) ’ “ 2(r^i) 


which exist and are distinct provided 

tG (1,9/5). (3.23) 

By direct calculations we get that t(j)'{x^'^) < 1 and t(/((xl^^) > 1, see Fig. [HI Hence, xl^^ is 
( 2 ) 

stable, whereas x* ^ is unstable. This means that 




■f / (2) 

if xi < x; , 


lim Xk 

/c—)-+oo 


„( 2 ) 


-|-oo 

\ 


■f (2) 

if xi = x; , 

■f ^ (2) 

if Xl > X* . 


(3.24) 


Note also that xi^^ —)• 1 and xi^^ —>■ oo as t ^ 1. 

Lemma 3.3. Assume that xi > x* , and hence x^ —>■ oo as A: —oo. Then the sequence 
{T{xk)}keN defined in fS. 1 6\) is not ergodic. 


Proof. The proof will be done by showing that the sequence of 


Ik = T23{xk) + T 32 (xfc), A; G N, 

is summable. This will imply that the sequence {T{xk)}k£N fails to have the infinite flow property 
in the sense of |28l Definition 2, page 1479]. This implies in turn, see |28l Theorem 1, page 1480] 
and also [25], the property in question. Standard linearization yields 

Xk+i — xl^^ > >r(xfc — xl^^), x := t())'(x®). 
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Figure 8 : Graphical solution of (13.211) 


which can be iterated to give 


Xk+l - > x^{xi - 


At the same time, by (I3.16p it follows that 

3 


Ik = 


Xfc + 3 


< 


Xfc — X: 


( 2 ) 


< 


_™( 2 )’ 


which is summable as x > 1 . 


□ 


By (I3.20p we obtain from the latter the following 

( 2 ) 

Corollary 3.4. For each xi > x\ , there exist e > 0 and the observable f such that, for each 
k gN, one finds I > k with the property d{Yi) > e. Hence, the model is in an ordered state. 

( 2 ) 

The unstable solution x* ' corresponds to the critical point which separates two basins of 
attraction, see (I3.24p . Note that the sequence {T{xk)}keN is still ergodic since all the entries of 
T(xl^^) = T(xfc), fc G N, are strictly positive. 

Now let us turn to the condition ()3.23l) . It yields 


P < f^{L) 


3 - x/5 

-s/5exp(L) — 3exp(—L) + 3 — 


As fi{L) is a decreasing function, the equation 'fi{L) = 1 has a unique solution 

L* = 7 ln^. (3.25) 

4 5 

For L < L*, one has 'fi{L) > 1, which means that t < 9/5 for all p G (0,1]. For such L and p, 
let KfiL,p) be the solution of the equation 

exp(4iF) = xl^^ > 1. 

Then ()3.24p can be rewritten in the form 



r 

vF :ijj 

if 

K < KfiL,p)- 


lim Xfc = < 

fc^ + OO 

^(2) 

if 

K = KfiL,p)- 

(3.26) 


, +00 

if 

K>AKfiL,p). 



For t = 9/5, we have xl^^ = xl^^ = 3, which corresponds to iF* = (ln3)/4. In this case, Xfc ^ 3 
if AT < iF*, and x^ —+oo if iF > iF*. 

For L > L* there exists p* = < 1 such, that for p G (0,p*), there exists KfiL,p) with 

the properties as in (|3.26p . For p G [p*, 1], the whole graph of tfi lies above the line tfi{x) = x, 
which means that x^ —>■ +oo for all initial xi > 0 . 

The results of the analysis just performed can be summarized in the following form. 
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Theorem 3.5. The Ising model on the graph based on Mi and described by the Hamiltonian 
h3.2i) has the following properties related to Definition \3.1\: 

(i) for L <{), it is in an unordered state for all values o/ K € M and p € [0,1]; 

(ii) for L G (0, L*] as in Hd.25\) and p € (0,1], there exists K^{L,p) > 0 such that the model 
is in an unordered state for K < K^{L,p), and in an ordered state for K > K^{L,p); for 
K = K^{L,p), the model is in the critical state; 

(Hi) for L > L*, there exists p* G (0,1) such that, for p < p*, there exists K,,{L,p) with the 
properties as in item (ii); for p G [p*, 1], the model is in an ordered state for all K. 

4 Concluding Remarks 

In this article, we introduce hierarchical random graphs based on motifs presented in Fig. [TJ The 
construction principles resemble those used in m- a nonrandom skeleton (hierarchical diamond 
lattice in m) is accompanied by random bonds. In our case they repeat the motif used in the 
construction. As a result, the motif appears at each hierarchical level. The construction is 
performed in a rigorous way and is illustrated by an informal description. The analysis of the 
node degree distribution in the constructed graphs is based on characteristic functions obtained 
in an explicit form. For p > 0, these functions are meromorphic for all motifs. This means that, 
for all our graphs, the node degree as a random variable has all moments with the property 
(n™) ~ C^ml. Thus, the degree distributions are intermediate as compared to the Poisson and 
scale-free cases. Such properties as clustering and small world property are studied only for 
p = 0,1. In particular, it turns out that for all motifs the small world property is absent for 
p = 0 and present for p = I. Thus, it would be interesting to hnd out how and at which value 
of p is emerges. In Theorem 13.51 we analyze phase transitions in the Ising model based on motif 
Ml. Unlike to m in our case the Ising model has no phase transition for p = 0, which manifests 
the difference between our construction and that used in m- We also show that, for L < L^, 
the model is in an unordered state whenever K = 0, i.e., the spin-spin interactions along the 
nonrandom bonds is absent. For L > and p > p*, the model is in an ordered state even for 
K = 0. We plan to study the phase diagram of this model in the {K, h)-plane in a separate work, 
where we also plan to consider such problems for the graphs based on the remaining motifs. 

Acknowledgment: This work was supported in part by the DFG through the SFB 701: ‘Spek- 
trale Strukturen und Topologische Methoden in der Mathematik’ and by the European Com¬ 
mission under the project STREVCOMS PIRSES-2013-612669. 


5 Appendix 

Here we give detailed calculations of the quantities from subsection [231 First we get the quantity 
in (12.91) : 






(5.1) 


l=l u=0 









X 



2 (g- 1)(1 -p)(l ^) 2p(fcg^ ^(l-g)+g-g^ ^) 







2 (g-l + (g-l)(fc-l)p) _ q^-\2q - 2 + 2p) - 2p 

qk-i j^i - qk-i I 
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Next, for the quantity in (|2.1ip . we have 

o A-i 


1=1 u=0 T \ / 


2(fc-l) 


+ E p+0 = 

y =0 ^ ^ 


4fc-i + 1 

= 4 + ifp-^±^ 
3 V 4^"^ + 1 

Now we calculate the quantity in (I2.13P 


k-l 3(1-1) 


(rik) = X] [ Z] 


i-l 




1=1 u =0 

4(1-1) 

+ Z (^+ 

u=0 

5(1-1) 

+ Z (^ + 

u=0 

2(k-l) 

+ Z (^ + 

i /=0 


1 + 4i-fc 


^)^p^(l-p)3a-l)-^ + 


4(/-l) 


1 + 41 -^ 

4 -' /5(/ - 1) 




'(l_p)4(/-l)-.' + 


1 + 41 -^ 

4 /2(A:- 1) 


3(fc-l) 

+ Z (3 + ^, 

i /=0 


4fc + 4 
2 

4^ + 4 


p^(l_p)5a-i)-^ 


+ 


fc-l 


+ Z(^ + 


u=0 


k-l 
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